Proof of Green's Formula 



Green's Formula: For the equation 



P{D)y = f{t), y{t) = for t <0 



(1) 



the solution for f > is given by 



y{t) = {f*w){t)= f f{T)w{t-T)dT, 



(2) 



where w{t) is the weight function (unit impulse response) for the system. 

Proof: The proof of Green's formula is surpisingly direct. We will use the 
linear time invariance of the system combined with superposition and the 
definition of the integral as a limit of Riemann sums. 

To avoid worrying about 0^ and f+ we will assume that f{t) is contin- 
uous. With appropriate care, the proof will work for an f{t) that has jump 
discontinuities or contains delta functions. 

As we saw in the session on Linear Operators in the last unit, linear 
time invariance means that 

y{t) solves P(D)y = f{t) =^ y{t - a) solves P(D)y = f{t - a). (3) 

Or, in the language of input-response, if y{t) is the response to input f{t) 
then y{t — a) is the response to input f{t — a). 

First we will partition time into intervals of width M. So, Iq = 0, ti = 
At, t2 = 2At, etc. 



Next we decompose the input signal f{t) into packets over each inter- 
val. The kth signal packet, fk{t) coincides with f{t) between tk and t^+i and 
is elsewhere 



= to ti t2 tk tk+1 

Figure 1: Division of the f-axis into small intervals. 
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Figure 2: The signal packet fk{t)- 
It is clear that f or f > we have f{t) is the sum of the packets 

m=Mt)+Mt) + ...+Mt) + ... 

A single packet /jt(^) is concentrated entirely in a small neighborhood 
of tfc so it is approximately an impulse with the same size as the area under 
fkit). The area under /^(f) !^ /(t^) At. Hence, 

fk{t)^{fitk)At)S{t-tk). 

The weight function w{t) is response to S{t). So, by linear time invariance 
the response to /^(f) is 

Vkit) ~ ifitk) AO w{t - tk). 

We want to find the response at a fixed time. Since t is already in use, 

we will let T be our fixed time and find i/(r). 

Since / is the sum of f^, superposition gives y is the sum of y^. That is, 
at time T 



y(T) = yo(T)+yi(T) 



(j{to)w{T - to) +f{ti)wiT -ti) + ..?j At 



(4) 



We can ignore all the terms where t^- > T. (Because then zv{T — t^) = 0, 
since T — t^ < 0.) If n is the last index where < T we have 

y(T) ^(^f{to)w{T - to) +f{ti)w{T - ii) + . . . + f{tn)w{T - t„)^ At 



t 
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This is a Riemarui sum and as At it goes to an integral 

y(T)= r f{t)w(T-t)dt 
Jo 

Except for the change in notation this is Green's formula (2). 

Note on Causality: Causality is the principle that the future does not affect 
the past. Green's theorem shows that the system (1) is causal. That is, y{t) 
only depends on the input up to time t. Real physical systems are causal. 

There are non-causal systems. For example, an audio compressor that 
gathers information after time t before deciding how to compress the signal 
at time t is non-causal. Another example is the system with input f{t) and 
output y{t) where y is the solution to y = f{t + 1). 
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